Skip to main content
Version: 2024.05

Charge system relaxation calculation

Relaxation for C-doped BN system

Input files

atom.config

    72
LATTICE
14.99639988 0.00000000 0.00000000
-7.49819994 12.98726326 0.00000000
0.00000000 0.00000000 20.00000000
POSITION
5 0.05556000 0.11111000 0.25000000 1 1 1
5 0.22221999 0.11111000 0.25000000 1 1 1
5 0.38888998 0.11111000 0.25000000 1 1 1
...
7 0.61111000 0.88889008 0.25000000 1 1 1
7 0.77777999 0.88889008 0.25000000 1 1 1
7 0.94444005 0.88889008 0.25000000 1 1 1
6 0.38888999 0.44444001 0.25000000 1 1 1

etot.input

 4  1
JOB = RELAX
NUM_ELECTRON = 288
IN.PSP1 = B.SG15.PBE.UPF
IN.PSP2 = C.SG15.PBE.UPF
IN.PSP3 = N.SG15.PBE.UPF
IN.ATOM = atom.config
RELAX_DETAIL = 1 100 0.01
Ecut = 50
Ecut2 = 200
MP_N123 = 1 1 1 0 0 0
E_ERROR = 0.0
tip
  1. NUM_ELECTRON is the total number of occupied valence electron in the system, one can use it to make the system charged.

B.SG15.PBE.UPF, C.SG15.PBE.UPF, N.SG15.PBE.UPF

tip

B.SG15.PEB.UPF, C.SG15.PBE.UPF, and N.SG15.PBE.UPF are the pseudopotential files.

Calculations

  1. You can submit PWmat tasks in different ways:
   mpirun -np 4 PWmat | tee output
Note

Run the command directly

   #!/bin/bash
#PBS -N SCF
#PBS -l nodes=1:ppn=4
#PBS -q batch
#PBS -l walltime=100:00:00

ulimit -s unlimited
cd $PBS_O_WORKDIR

mpirun -np 4 PWmat | tee output
Note

Submit the task with a pbs script

  1. For RELAXATION calculation, files to watch during relaxations.

RELAXSTEPS: concisely reports about the atomic relaxation steps. A typical RELAXSTEPS file looks like:

   It=   -1 CORR E=  -0.1267760205584E+05  Av_F=  0.61E+00 M_F=  0.41E+01 dE=  0.48E-03 dRho=  0.40E-04 SCF=    33 dL=  0.00E+00 p*F=  0.00E+00 p*F0=  0.00E+00 Fch=  0.00E+00
It= 0 NEW E= -0.1267351064854E+05 Av_F= 0.20E+01 M_F= 0.12E+02 dE= 0.21E-03 dRho= 0.38E-04 SCF= 18 dL= -0.37E+00 p*F= 0.29E+02 p*F0= -0.90E+01 Fch= 0.93E+00
It= 1 CORR E= -0.1267810094490E+05 Av_F= 0.17E+00 M_F= 0.71E+00 dE= 0.56E-04 dRho= 0.44E-04 SCF= 11 dL= -0.95E-01 p*F= 0.12E+00 p*F0= -0.90E+01 Fch= 0.10E+01
It= 2 NEW E= -0.1267804036855E+05 Av_F= 0.31E+00 M_F= 0.17E+01 dE= 0.26E-03 dRho= 0.42E-04 SCF= 12 dL= -0.83E-01 p*F= 0.36E+01 p*F0= -0.24E+01 Fch= 0.92E+00
...
It= 16 CORR E= -0.1267817922913E+05 Av_F= 0.21E-02 M_F= 0.67E-02 dE= 0.54E-04 dRho= 0.14E-04 SCF= 2 dL= -0.55E-03 p*F= -0.41E-02 p*F0= -0.60E-01 Fch= 0.72E+00
It= 17 *END E= -0.1267817922913E+05 Av_F= 0.21E-02 M_F= 0.67E-02 dE= 0.54E-04 dRho= 0.14E-04 SCF= 2 dL= -0.55E-03 p*F= -0.41E-02 p*F0= -0.60E-01 Fch= 0.72E+00

REPORT: more detailed information about every electronic and ionic step

  1. Other important files.

final.config: holds the structure of the last ionic step, the structural result (also very important for restarting a relaxation)

MOVEMENT: holds the structures of every ionic step during relaxations.